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ABSTRACT 



The reorder point and reorder quantity for a multi- 
echelon inventory system consisting of two levels were 
determined through the use of a mathematical model and a 
computer simulation was used to verify the results. 

The main echelon supported two lower echelon stock 
points and reordered using a continuous review inventory 
policy. The two lower echelon stock points operated under 
periodic review policies. 

The measure of effectiveness used was to minimize total 
system costs subject to a constraint on the maximum number 
of back orders per year. 

The results from the mathematical model were used as 
input values for the simulation and measures of effective- 
ness were compared. An alternative procedure was proposed 
and simulated; and the results of the three products were 
compared. 
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I. 



INTRODUCTION 



Despite the fact that all military supply systems and 
many large civilian corporations have multi-echelon and 
multi-location inventory systems, few multi-echelon inven- 
tory models are currently being used. The policies that 
have been used are the single echelon and single location 
inventory policies that attempt to minimize the total 
variable cost of a single location ignoring total system 
cost. The difficulties of optimizing a multi-echelon and 
multi-location inventory system are most likely due to the 
complexity of demand and the interdependency among the units 
in different echelons. 

A mathematical model for a multi-echelon and multi- 
location inventory system is developed in Chapter II. The 
objective of this model is to minimize the total annual 
variable costs subject to a maximum allowable number of 
back orders per year for the entire system. This model 
consists of one first echelon location operating under a 
continuous review policy and two second echelon systems 
which operate under periodic review policies. An example 
is presented to obtain numerical results. 

In Chapter III, a computer model is developed to simu- 
late the multi-echelon, multi-location system. This simu- 
lation model was used to check the results obtained from 
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the mathematical model and to demonstrate the interaction 



between echelons. Graphical plots are generated which show 
the inventory position of each entity in the model. 

In Chapter IV, comparisons are made of the simulation 
results and the results obtained from the mathematical 
model. The comparisons suggest some modifications to the 
operating policy. 

The system was simulated again using the modified 
operating policy and substantial improvements in the 
effectiveness of the multi-echelon system were observed. 
These results are described in Chapter V. 

Chapter VI summarizes the results of the research and 
concludes with some suggestions for additional research. 
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II. A PROBABILISTIC MATHEMATICAL MODEL FOR MULTI-ECHELON 
INVENTORY SYSTEM 



A. DESCRIPTION OF MODEL 

In this model there are three systems. One of these 
systems represents the highest echelon and is called the 
Main System. The other two, at separate locations, are 
the lower echelon and are called System One and System Two, 
respectively. Each system carries its own inventory, and 
receives random demands. System One and System Two are 
dependent upon the Main System but independent of each 
other. System One and System Two can only be resupplied 
by the Main System. In other words, they cannot order 
from any other suppliers. The Main System replenishes 
stocks by placing orders to external suppliers. 

B. OBJECTIVE OF MODEL 

The objective of this model is to minimize the total 
expected annual variable costs of three systems subject 
to a specified expected number of total back orders per 
year. In fact, this is the same as minimization of total 
yearly variable costs subject to a specified minimum level 
of customer satisfaction. 

Based on these objectives, decision variables for each 
system will be calculated to achieve optimum levels for 
the entire multi-echelon supply system. 
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C. ASSUMPTIONS OF THE PROBABILISTIC MATHEMATICAL MODEL 



1 . Main System 

It is assumed that the Main System under considera- 
tion consists of. a single installation which utilizes trans- 
action reporting. This system is governed by a (Q, r) type 
policy with back orders. 

The assumptions in addition to the continuous review 
assumption are: 

a. The cost of operating the information processing 
system is independent of Q (reorder quantity) and r (reorder 
point) . 

b. The unit cost C of the item is a constant in- 
dependent of Q. 

c. The back-order cost is constant (n) , per unit 
back ordered regardless of the length of time the back- 
orders exist. 

d. There is never more than a single order out- 
standing. This assumption implies that when the reorder 
point is reached, there are no orders outstanding; there- 
fore, the inventory position is equal to the net inventory. 
Thus, the reorder point will be the same regardless of 
whether it is based on the inventory position or net 
inventory. 

e. Procurement lead times are independent and 
identically distributed (i.i.d) random variables with a 
gamma distribution. 
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f. All variables are treated as continuous. 

g. The demands are Poisson distributed with the 
mean number of demands per year a constant A . 

h. The reorder point, r, which is based on the 
inventory position is positive. 

With a back orders constraint, it is infeasible to 
wait until back orders exist before placing an order. 

Because of this and assumption (d) there will be no back 
orders outstanding at the reorder point. As was discussed 
in assumption (d) , at the reorder point, the inventory 
position is equal to the on-hand inventory. 

For this model, any one of the three inventory 
levels on hand, net, or inventory position can be used to 
define the reorder point; and the reorder point has the 
same value for any one of them. It should also be noted 
that to use the on-hand level it must be assumed that 
after an order arrives, it is sufficient to fill all back 
orders and raise the on-hand inventory level above the 
reorder point. If this ever failed to happen, the reorder 
point would never be reached again and the system would 
continue to accumulate back orders indefinitely. 

When the reorder point is thought of in terms of 
the inventory position of the system, then assumption (d) 
guarantees that the on-hand inventory will always exceed 
the reorder point when an order arrives; otherwise, it 
would not be possible to have only a single order outstanding. 
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2 . System One 



This system, which is at the lower echelon, is 
governed by a periodic review policy. The operating 
doctrine is the most widely used type of periodic review 
which is the order up to R policy. All demands which 
occur when the system is out of stock are back ordered. 

For this periodic review system the time between 
reviews will be denoted by T, and at each review time a 
sufficient quantity is ordered to bring the inventory 
position of the system up to a level, R, regardless of the 
amount of on-hand inventory. This policy dictates that 
at review times even if the inventory position is R-l, 
only one item must be ordered to bring the inventory 
position up to R, ignoring the high cost of placing the 
order. Despite its appearance of being illogical, this 
assumption simplifies the formulation of the system and 
is very unlikely to occur on high demand items. When an 
item experiences zero demand in a review period, there is 
no need for order because the inventory position is 
already at R. 

The other assumptions are: 

a. The cost, J, of making a review is independent 
of the variables R and T. 

b. The unit cost, C, of the item is constant and 
independent of. the quantity ordered. 
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c. Back orders are incurred only in very small 
quantities. This implies that when an order arrives, it 
is almost always sufficient to meet any outstanding back 
orders . 

d. The backorder cost is constant, IT, per unit 
back ordered regardless of the length of time the back 
order exists. 

e. Procurement lead times are i.i.d random 
variables with a gamma distribution. 

f. Orders are received in the same sequence in 
which they were placed. It should be noted that for (Q, r) 
models, the two assumptions that orders were received in 
the sequence placed and that lead times are i.i.d random 
variables could not both hold rigorously, since there 
exists a positive probability that two successive orders 
could be separated by an arbitrarily short time interval. 

In this model, orders can never be more closely spaced 
than by an interval of length T. If T is large enough, 

it is possible, provided that there is a sufficiently 
small range of variation in the lead time, that both 
assumptions hold simultaneously. 

g. The demands are Poisson distributed with 

mean, A. 

h. All variables are treated as continuous. 
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3 . System Two 



This system is identical to, but independent of. 
System One. 

D. ANALYTICAL SOLUTION TO THE MODEL 

As was mentioned previously, the objective of the model 
is to minimize total system costs subject to a constraint 
on the number of back orders per year. 

The determination of various annual costs of each 
system may be made independently and placed into a common 
cost formula. Instead of doing this independently, it is 
better to find a main system cost expression. System One 
cost expression and System Two cost expression, and then 
to add them to each other to determine the total cost 
formulation and then use this formulation as the objective 
function of total system. The systems are actually tied 
together through the constraint on the number of back 
orders . 

For the following notation, the subscript m will 
indicate that the variables being subscripted belong to 
the Main Sytem where 1 and 2 indicate System One and 
System Two, respectively. 

The costs that are of interest in this model for each 
system are the cost of placing an order, the cost of 
carrying inventory and back order costs. 
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1 . Main System 



As was discussed previously, the Main System has 
a transactions reporting policy or (Q, r) model. 

In a continuous review system, a period is defined 
as the length of time between the receipt of two successive 
procurements . T his time peri od is a ran dom variable. 

Because the procurement lead times are random variables, 
the number of demands for a fixed time are random variables; 
and the number of items demanded per demand are also random 
variables. The review period is also a random variable. 

In the following material, we itemize the costs. 

a. Procurement Cost 

= Expected number of demands per year 

Q m = Order quantity 

A 

-r^ = Average number of procurements per year 

^m 

A = Procurement cost per cycle 
m 

A 

rp- A = Procurement cost per year 

y m 

b. Inventory Carrying Costs 

Because of randomness, it is possible to accumu- 
late a large number of back orders at the end of a cycle. 

To prevent this from occurring, it is advisable to provide 
a safety stock to buffer the system from excessive numbers 
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of back orders. The safety stock is the expected amount 
of stock on hand when an order arrives. The actual amount 
of stock on hand when a shipment arrives is clearly random. 

S = Mean value of on-hand stock when an order 
m 

arrives 

After an order arrives, the expected on-hand inventory in- 
creases to Q+S and is reduced to a value of S on the average 
just before the next order arrives. Therefore, the average 
on-hand inventory per cycle is: 



(Q +S ) 
m m 



+ -f = 8 + 

2 2m 



To write S in terms of the reorder point, r , 
m m 

let us first assume that the lead time x is fixed. Let 

m 

E (x;r ) = r - x be the net inventory 
x m m 

at the time an order arrives, 

where x 

is the number of units demanded in lead time x . Then, 



00 

S = [Net inventory] = / E (x;r ) f(x;x ) dx 

m t o Tm m m 

m 



where f(x;x ) = Density function of demand in time t 
m d m 



oo CO 

S m = / (r -x) f ( x ; x ) dx = r / f(x;x ) dx 
m o m m m o m 

00 

- . / f (x-;T ) dx 
o ' m 
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S m “ r m _tJ m where y m = Expected lead time demand. 

Then ^ + S = ^ + r - y 
2 m2 m m 



I = Inventory carrying charge 
C = Cost of an item 



Total holding cost/year = IC (r - y + 9^) 

u mm2' 



Stockout Costs 



0 if x-r <0 

Let us define n (x;r ) = { m 

X m « fm ‘ 

m x-r if x-r ^0 

m m 



where n T (x ; r m ) is the number of back orders per cycle, 
m 

If n m (r m ) = expected number of back orders per cycle, then 
n m (r ) = /°° n_ (x* r ) h (x) dx = /°° (x-r ) h(x) dx 

mm o m o m 

= f°° xh(x) dx-r m H(r m ) where H(r m ) = P [X>r m ] 
r 

m 

and h(x) = marginal distribution of leadtime demand. 
Therefore, the expected number of back orders/year = 



f l r xh(x) dx - r m H(r m ) ] 
v m m 



and the expected cost of back orders/year = 



n A 

-Sp if xh(x) 
y m m 



dx-r H ( r ) ] 
m m 1 
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All the terms in the average annual variable 



cost K have now been found: 
m 



K 



m 



A 

10 



+ IC 




r 



m 



iu_ 

y m ] + Q 



m 



m 



[ /°° xh (x) dx-r H (r ) ] 
r mm 

m 



2 . System One and System Two 

It was stated in the assumptions that System One and 
also System Two both follow a periodic review policy with an 
order up to R stockage policy. Since System One and System 
Two are identical to each other; the equations will be 
derived only for System One. 

For convenience, a period is assumed to be the time 
between the receipt of two successive orders rather than 
between the placement of two successive orders. Costs are 
described as follows: 

a. Ordering and Reviewing Costs 
J = Reviewing cost/cycle 

A = Ordering cost/cycle 

Jl+Ai 

Ordering and reviewing cost/year = — = 

L 1 

where Tj_ is the period length defined in units, of years. 

b. Inventory Carrying Cost 

The expected net inventory just prior to the 
arrival of an order is R^ - - X^T^, where = mean demand 

during lead time. 

The mean rate of demand remains constant over 
time and the expected demand per period must be the expected 
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If the expected net inventory 



amount ordered, i.e., A^T . 
immediately after the arrival of a procurement is R^-y^, 
it is therefore ~ V 1 ~ just prior to the arrival 

of a procurement. 

The expected unit years of storage incurred 
per period is 

A T 

T 1 (R i - V + 7 < R i - Hi - X lTl )] = Tl [ Rl - Vl - Jj-i] 

AfTi 

and average inventory carrying cost/year = IC [R^ - y 1 - — ^ — ] . 
c. Stockout Costs 

First we assume the case where the procurement 
lead time is constant t. An order placed at time t will 
arrive in the system at time t + x , and the next procurement 
will arrive in the system at time t + x + T-^. After the 
order is placed at time t, the inventory position of the 
system is R^. It is necessary to compute the expected 
number of back orders occurring between t + x and t + x + T^. 

A back order will occur in this period under assumption c if 
and only if the demand in the time period x + T. exceeds R . 
Assumption c also assures that after the arrival of the order 
placed at time t, there will be no remaining back orders, 
and therefore they must all occur between times t + x and 
t + x + Tj. Consequently the expected number of back orders 
incurred per period is 
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R / (x-R^) f(x; x + T^) dx where 

f (x; t + T-^) = Demand distribution during time x + T^. 

When lead time is random with density g (x^) 
with x • „ and x being lower and upper limits respectively 
andx^andx 2 , the lead times for the orders placed at times 
t and t + , respectively, the expected number of back 

orders incurred per period is: 

T T 

max max 

oo 

T / T / R / (x-R^) f (x ; t 2 +T 1 ) g (x 2 ) g (Tj^) dx dx 2 dx 1 
‘min min 1 

A 

= r /°°(x-R 1 ) h (x; T^) dx 

/N T 

where h(x; T^) = I maX f(x; t g ( t dx ^ which is the 

min 

demand distribution during time x 2 + when lead time is 
a random variable with density function g(x 0 ). The average 
number of back orders incurred per year is; 

E l (R l , T l ) = r^"°° (X-R 1 ) h(X; T 1 ) dx and the average 

back order cost per year equals toII^E^ (R 1 ,,T 1 ) and 

ri^Cr) = T ]_ E ]_ ( R ^ »T 2 ) is the expected number of back orders 

per period. 
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Finally, the annual variable cost of System One 

is : 

Li X..T.. 

= — + IC [R x - y 1 - + n x E x (R^T.^ and 

likewise the annual variable cost of System Two is: 

K 2 = Tj + IC [ R 2 y 2 ~ 2 ^ + n 2 E 2 ^ R 2' T 2* 

where Li = J . + A . . 

1 i 

3 . Objective Function and Constraints 

The objective function of the model consists of the 
total annual variable costs of each system. Therefore; 



K = K + K. + K~ which is equal to the minimization 
m 1 2 ^ 



of; 



x q nx 

K = pj— ■ A + IC [-£ + r - y ] + n (r) + 1 

Q m 2 mm m 



m 



m 



h 

T, 



+ IC [R x - M 1 - ~Y±] + ^ n 1 (r) + ^ + IC [R 2 -y 2 — — J 
n 2 - 

+ m=- n 2 (r) 

X 2 



subject to n (r) + p- (r) + p- n 2 (r) <b 
y m 1 2 



where b is the specified total maximum number of back orders 
per year for the entire system. 

The problem at hand is to calculate the optimum 
values of Q m , r m , R 1 , Tp R 2 , T 2 . Since the objective 
function and the constraint are non-linear functions of the 
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decision variables, we solve the problem using the Lagrange 
multiplier approach. 

4. Optimum Values of Operating Variables 

After including the Lagrange multiplier in the 
formulation, the new objective function becomes: 

Minimize 



L = pr- A+IC[-J^ + r - u ] + 

m 2 m m J Q 
m ^ 



- . . , h 

„ R m r + t7 
m 1 



ApTp-i 

+ IC f R l - H - -Vl + T - V r > + ^ + ic[r 2 - y 2 - 

l ^ 



+ T_ ^2 (r) 9 ^ ( Q ^m (r) + T, ^l (r) + T- ^2 (r) J 

2 m 1 2 



- b] 



where 0 is the Lagrange multiplier. It is through this 
Lagrange multiplier that the three systems are linked together 
mathematically . 

The optimum values of the unknown variables can be 
found by taking the derivatives of the objective function 
with respect to Q^, R^, R^, T^, R£ , , 9; equating them to 

zero; solving the equations simultaneously; and ensuring the 
Kuhn-Tucker conditions are satisfied. 

The derivatives were taken with respect to Q , r , 

R^ , R£ , 0. A different procedure was utilized to find T^ 
and T 2 • 
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The derivatives: 



_ 3L 

J Q 3Q 
m 



A A 

m m , IC 
m 

A 



~Y~ n_(r) (n 
m m 

nn 



- 0) = 0 



m 



lr = IC ' of V r > <"m - 9 > - 0 

m 



L R x = - IC - ^ < n i - S) - 0 



3L 1 

L R 2 = 3Rj Ic " T^ H 2 (R 2' T 2 ) (n 2 9) = 0 

L e = H = t ( q^ (r) + F" V r) + T~ ^2 (r)) - b ] = 0 

m 1 2 



The Kuhn-Tucker conditions : 



L ^0 Q . L = 0 

Qm m Q 

m 



Q m * o g £ o 



L 5i 0 r L =0 

r m r 

in m 



r > 0 0 . g = 0 

m * 



\ 4 0 R i- L Rl = 0 



R 1 ^ 0 0 < 0 



L i 0 
R 2 



R 2‘ L R 2 = 0 



R 2 * ° 



where g = gE n m (r) + fj- n^r) + n 2 <r) - b 



Therefore the following Kuhn-Tucker conditions must 
be satisfied: 
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A m A m IC A m - (r) (n - 

— qT" + T ‘ Q" n m 
m m 



ic - ^ h (r) (n - e) a o 

Q m m 

m 



ic - ~ h 1 (r 1 , t x ) (n 1 - 0: 



IC " — 2 H 2 ( V V (II 2- ( 



X A IC Q X 
m m nn m - 

% 2 Q m m 



X r 

IC r - H (r) (II - 9 

m Q m m 

m 



ic r i - ^ s i (E i- t i> <n i 



IC R 2 " T 2 ®2 * R 2' T 2* (K 2 



5 0 



r _ 5 0 
m 



R x * 0 



R 2 > 0 



^m (r) + T, ^l (r) + T 2 ^ 
1 2 



0 [<f V r > + ¥: V r > + ?; 

m i* ^ 

9^0 



9) * 0 

I } 0 

9) :> 0 

(n - 0) = 0 

m 

) = 0 
- 0) = 0 
- 6 ) = 0 



2 (r ) - b < 0 
• n 2 (r) - b] = 0 
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From these equations and conditions the optimum 
values of the unknown variables can be simplified to the 
following form and the exact values can be found by 
solving these equations simultaneously: 

IC Q 

Hm(r) = 

q = / 2 \ [a r~n Hi _ 0) i 

m m L m + m m J 

IC 

IC T, 

H 1 (R 1' T l> = tv^t 
IC T 2 

H 2 (R 2' T 2 ) “ (n 2 -0) 

The procedure to find and T 2 and subsequently 
and R 2 is through iteration and trial and error. This is 
not unreasonable. In realistic cases other considerations 
not modelled here usually dictate the length of a period. 

In most cases the systems control a large number of different 
items and the same period length is used for each item. 

Thus, the period length is usually some convenient calendar 
or financial period, such as a month or a quarter. The 
procedure we use is shown in the following example. 

E. EXAMPLE PROBLEM 

1. Additional Assumptions 

For this example it is assumed that each system has 
Poisson arrivals independent of each other. The customers 
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demand only one item at a time. System One and System Two 
are resupplied only from the main system and the main system 
is resupplied from outside suppliers. To be more realistic 
it is also assumed that the procurement lead* times have a 
gamma distribution with different mean and variances. 

The X's will be assumed daily demand or arrival 
rate per day rather than yearly values. 

2 . Input Values 

Before providing the input values, some clarifications 
must be provided regarding the demand at the main system. 

When it is stated that = 4/day, X ^ = 3/day and = 1/day. 
These represent the mean number of demands which arrive each 
day directly at the respective systems. However, since all 
demands at Systems One and Two eventually must filter up to 
the main system, the cumulative demand at the main system 
has an expected value of + X^ + ^ m = The ma i- n system 

supports not only Systems One and Two, but also has its own 
customer demands. The model assumes that the demand at the 
main system is the superposition of the direct demands at 
the main system and Systems One and Two. However, since the 
actual demands placed on the main system in the lower echelon 
are batched, the demand variability is much greater than 
would be expected by the superposition process. We will 
evaluate the seriousness of our assumption about the demand 
process at the main system with the simulation model 
described in the next chapter. 
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In this example the input values for each system 



are ; 

a. Main System 

X m = 8/day (considering the individual customer 
arrivals and the remainder being Systems 
One and Two average demands) 



c = 


$50/item 


I = 


0.23 


n m = 


$5 


A m = 


$500 



b. System One 



h ■ 


3/day 


C = 


$50/item 


I = 


0.23 


n i = 


$5 


A 1 " 


$100 


J 1 = 


$100 



c. System Two 



X 2 


1/day 


C = 


$50/item 


I = 


0.23 


n 2 " 


$5 


A 2 = 


$100 


J 2 = 


$100 
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The total system is expected not to exceed 



b = 175 total number of back orders at the end of a year. 

3 . Solution Procedure and Results 

It was assumed that for each system, procurement 
lead times are gamma distributed random variables. 

The gamma distribution is defined as: 

— 

g(t) = t a ^ e ^ with 

r(a) s a 

y = a • 8 and 
Var = a. 6^ 

\ 

We need the lead-time demand distribution for the 
case in which lead times are gamma distributed and the 
process generating demands is Poisson; This is derived 
below. 

The Poisson distribution is 
, , , > x -At 

f (x) = — — — and the demand distribution 

x J 

during lead time is: 



f(x;x) = / f(x) g(t) dt 
o 



= / 
o 



Ut) x e - Xt t 0 ' 1 e' t/B 



X 



r(a) e a 



dt 



* X , /" t- 1 + x e - ( 6 + X)t dt 



r (a) B a x 1 0 
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, X 



r (a) ea x * 



p (a+x) 

<F +X) 



,1 ,,a+x a-l+x (g+A)t 
r co ( f +A ) t e 

a+x o r ( a .+ x ) dt 



,» <| +A)“ +X t“- 1+x e" <e +X>t 

0 nS+J) * = 1 



SO 



f ( x ; T ) = 



x 



T (a+x) 



r (a) g 1 ^ xl ,1 ,.a+x 

( g +X) 



- T (a+x) 



r(a) x : l ( i +A) 



] x ] a 
l +x 



a+x-1 , -r 

< , ) ] x Et-S— ] a 

b x t+ x 



a+x-1 

X 



= ( ) (l-p) x p a 



_ a (1-p) 
M P 



a ( 1- 



1/6 



1/g + X 
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1/6 + A 



lag 



Var = ^ P J 
2 

P 



aa -T7ek ) 

1/6 

1/6+A 



Aag (1+Ag) 



In fact this is a Negative Binomial distribution 



with 



p 1/g+A 



The 


parameters 


a and 8 


for each 


system 


are : 


a = 
m 


: 12.8 


a i 


= 5.43 


a 2 


= 2.5 


S m = 


: 3.125 




= 4.6 


S 2 


= 6. 



The main system has a Negative Binomial distributed 
lead time demand with parameters: 

= 8- ( 12.8) (3.125) = 320 

var m = 8- (12.8) (3. 125) [1+8 • ( 3 . 125) ] = 8320 

System One has the parameters: 

V 1 = 3- (5.43) • (4.6) = 74.934 

var x = 3- (5. 43) (4. 6) [1+3(4. 6)] = 1109.0232 

System Two has the parameters: 

U 2 = 1* (2.5) -6 = 15 

var 2 = 1 * (2.5) *6 [1+1(6)] = 105 

Because the negative binomial is computationally 
intractable, we use the normal approximation for the calcu- 
lation of lead time demand probabilities. The normal 
distributions are assumed to have the same mean and variance 
as the negative binomial distributions they replace. 

To solve this example problem, it is necessary to 
make an initial estimate for the Lagrange multiplier 6 
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which must be less than or equal to zero. After the initial 
estimate if the total number of back orders per year exceeds 
b, then the absolute value of 0 should be increased gradually 
until the number of back orders converges to b. 

For this problem 9 = -1.5 works very well. For the 
main system the formulas are: 






H m (r) = (TT^TX " P [X S r m! 



m 



m 



0, is calculated to be 550 and 
in 

H (r) = P fX>r 1 = 5S5Q (0.23) * 50 q 3332455216 
m' r) L *" r m J (5+1.5) 2920 u • 444243:>2±b 

From the inverse standard normal distribution TI-59 
calculator program (Appendix G) , the reorder level at the 
main system is found to be: 



r = a .(0.4305333395) + y = 359.2706827 
mm m 



V r > - <V r m> V r > + a m 

m 



where 0 is the functional value of standard normal distribu- 
tion at 0.4305333395 

n (r) = 20.08139334. Using this value in the Q 
m m 

formula yields 
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568.84 



Q m = / 2T292Q [500 + 20.08139334 (5 + 175) 

(0. 23) -50 “ 

which is not equal to the first estimated value. If Q is 

m 

taken 569 then: 



V r > 



569 (0.23). 50 
(5 + 1.5) 2920 



0.3447576396 



r m = a m ( 0 * 399071201) + 0 m = 356.4008941 

n_ (r) = 21.05458646 
m 

Using this value in the Q formula yields 

m * 

Q = /2.2920 [500 + 21.05438646 (5 + 1.5) = 568.69 

(0.23) -50 



which is very close to the initial value. A summary of 
the results for the main system follows: 

Q = 569 

m 

S m = 36.400891 

m 

n m (r) = 21 . 05438646/period 

Fi m (r) = 108.047115/year 

K = $6796 . 50/year 

m 

To solve for the optimum values for Systems One and 
Two a different procedure is followed. The total annual cost 
of Systems One and Two is a convex function of the period length 
T. For different values of T there are different total annual 
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cost values. The minimum of these values is the optimum 
total annual variable cost and the corresponding T and R 
values are the optimum operating values. A TI-59 program 
was written to perform the line search for the best value 
of T. The program is found in Appendix G. User information 
and the features of the program are also in the same appendix. 
This program evaluates the R value, safety stock, back 
orders per period, back-orders per year, annual reviewing 
and ordering cost, annual holding cost, annual back order 
cost, and finally, the total annual cost. 

Using this program the optimum values for Systems 
One and Two are: 



R 



1 



T 



1 



h-L (r) 

r\ 1 (r) 



K 



1 



= 306.809389 

= 2.39 months (72.6958 days) 
= 13.78788901 
= 8.670289086/period 
= 43.53283223 
= $2634 . 41/year 



R 2 = 134.9807003 

T 2 = 4.07 months (123.7958333 days) 

S, = 0 

ri 2 (2) = 8.13480099/period 



34 



n 2 (r) = 23.98467122/year 

K 2 = $1377. 56/year 

The total back orders per year are: 

n m (r) + (r) + n 2 (r) = 108.047 + 43.532 + 23.984 

= 175.56 

which is the same as b = 175.56. 

The total annual variable cost is: 

K = K + K, + K_ = 6796.50 + 2634.41 + 1377.56 
m 1 2 

K = $10808.47. 
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III. COMPUTER SIMULATION FOR MULTI-ECHELON MULT I -LOCATION 
SINGLE ITEM INVENTORY SYSTEM 



To check the analytical results, a computer simulation 
was written which uses the same operating assumptions and 
input parameters that were made for the analytical model. 
This model, however, simulates the real world more accu- 
rately since some of the simplifying assumptions required 
to obtain analytical results were not necessary in the 
simulation. Also, the demands placed at the main system 
from the lower echelon systems were batched as in the real 
world. 

A. DESCRIPTION OF THE SIMULATION MODEL 

As in the analytical case, there are three systems in 
the simulation model; the main system. System One and System 
Two. The flow charts of this program are in Appendix A. 

1 . Main System 

This system uses a continuous review policy. When 
the stock level reaches the reorder point, it orders the 
amount Q. The decision variable for reordering is the in- 
ventory position. When an order is placed, the inventory 
position increases by an amount of Q. If the new inventory 
position is less than the reorder point, the system places 
an additional order for another amount of size Q. Then the 
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inventory position increases with one more Q. This con- 
tinues until the inventory position exceeds r. The order 
policy is thus (nQ,r) . N is the smallest integer that will 
make the inventory position higher than the reorder point. 

There are independent customer arrivals to the 
main system and also group demands that are placed by 
System One and System Two when those systems replenished 
their stocks at the end of their periods. 

The number of units demanded per requisition is 
one. The number of units demanded for resupply to the lower 
echelon systems is random depending on the demands at the 
lower echelon systems and the parameters R^, R2. The 
program is capable of allowing geometrically distributed 
quantities demanded per requisition. Back orders at the main 
system are satisfied by filling the back orders to individual 
customers first followed by filling any back orders due to 
System One and System Two on a first-come, first-served 
basis. 

When a demand occurs from the lower echelon systems 
for which there is not sufficient on-hand inventory at the 
main system, the maximum amount is filled and the rest is 
put into the back order queue. As soon as a shipment 
arrives at the main system, the back orders are filled. 

The times between the customer arrivals are independent 

of each other and exponentially distributed. The lead 

times are also independent of each other and gamma distributed. 
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2 . System One and System Two 

Both systems operate identically. The only differ- 
ence might be in the values of the system parameters. They 
have the same periodic review policy which is: order up to 

R at each review time. The program is also capable of making 
a decision at each review time regarding the placement of 
an order if the inventory position is less than or equal to 
a threshold value, is an (r,R) policy. If r is taken to be 
R-l, then the system orders up to R at each review time even 
if there is only one demand in a period. 

The times between customer arrivals are independent 
of each other and exponentially distributed. The number of 
units demanded per requisition is one. The lead times are 
independent of each other and gamma distributed. 

B. HOW TO USE THE PROGRAM 

Two simulations are given in Appendix C and Appendix D. 
The programs differ primarily in the type of output that 
is generated. The program given in Appendix C produces a 
Versatec plot output showing the inventory position for 
the main system and the lower echelon systems. A simulated 
period of up to 4 or 5 years can be run with a Class K job. 

On the Versatec output for each system, there are two 
plots: one plot has a small triangle at each point repre- 

senting the inventory position. The second plot represents 
the net inventory. The first figure shows the main system 
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values, the second figure shows the System One values and 
the. third figure shows System Two values all plotted versus 
time. 

The second program can be used to simulate operations 
over arbitrarily long periods of time. There are no dimen- 
sion restrictions. This program gives the net inventory 
at the end of each period for all the systems; total demand 
in the periods of both System One and System Two; demand 
during a lead time for the main system; average on-hand 
inventory; and average number of items short per day for 
each system. 

The input required by the programs is described by the 
variable definition list given in Appendix H. The programs 
use the IMSL subroutines for random number generation. 

C. RUNS 

1 . Starting Conditions 

There were two different starting conditions entered 
for the main system. In the first case the inventory posi- 
tion and net inventory were set initially to equal the order 
quantity. In the second case, the inventory position and 
the net inventory were set equal to the reorder point plus 
the order quantity. For Systems One and Two the inventory 
positions and net inventories were set initially at . 
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2. Results of Runs 



In all runs the length of time was 10 years. The 
average yearly results were obtained by dividing the result 
for the 10 year period by 10. 

Using the policy parameters determined by the mathe 
matical model and the first set of starting conditions the 
following simulated results were obtained, 
a. Main System 

The number of orders in 10 years = 51. 

The number of back orders in 10 years = 3026. 

The average on-hand inventory over 10 years = 

313.026. 



The average safety stock when an order arrives 
56.54 units. 

The average annual costs are: 

Ordering costs: $2,550.00 

Holding costs: 3,599.80 

Back Order Costs: 1,513.00 

Total Cost $7,662.80 



b. System One 

Number of back orders in 10 years = 587. 

The average on-hand inventory over 10 years = 
114.482 units. 
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The average annual costs are: 

Ordering and Reviewing Costs: $1,004.18 

Holding Costs: 1,316.54 

Back Orders Costs: 293. 50 

Total Costs $2,614.22 

c. System Two 

Number of back orders in 10 years = 287. 

The average on-hand inventory over 10 years = 

57.618 units. 

The average annual costs are: 

Ordering and Reviewing Costs: $ 589.68 

Holding Costs: 662.61 

Back Orders Cost: 143.50 

Total Costs $1,395.79 

(1) Total results 

The total number of back orders for the 
whole system per year is 302.6 + 58.7 + 28.7 = 390. 

The total annual variable cost for the 
whole system is $7,662.80 + $2,614.22 + $1,395.79 = $11,672.81. 

The model was run again with the second set 
of starting conditions and the same set of policy parameter 
values. The simulation results are summarized below: 

(a) Main System 

The number of orders in 10 years = 50. 

The number of back orders in 10 years = 

1665. 
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3,989.78 

832.50 

$7,322.28 



The average on-hand inventory over 
10 years = 346.937. 

The average safety stock when an order 
arrives = 77.46. 

The average annual costs are; 

Ordering Costs: $2,500.00 

Holding Costs: 

Back Order Costs: 

Total Costs: 

(b) System One 

Number of back orders in 10 years = 582. 

The average on-hand inventory over 
10 years = 118.696 units. 

The average annual costs are; 

Ordering and Reviewing Costs: $1,004.18 
Holding Costs: 1,365.00 

Back Orders Costs: 291.00 

Total Costs: $2,660.18 

(c) System Two 

Number of back orders in 10 years: 314. 

The average on-hand inventory over 10 years 
57.401 units. 

The average annual costs are; 
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Ordering and Reviewing Costs: $ 589.68 



Holding Costs: 



660.11 



Back Orders Costs: 



157.00 



Total Costs 



$1,406.79 



(2) Total Results 

The total number of back orders for whole system 
per year is 166.5 + 58.2 + 31.4 = 256. 

The total annual variable cost for whole system 
is $7,322.28 + $2,660.18 + $1,406.79 = $11,389.25. 

3 . Comparison of Both Starting Conditions 

It is obvious that if both results are compared there is 
a considerable decrement on annual back orders but not a great 
difference in annual variable costs. 



390-256 

390 



0.3436 



11672.81 - 11389.25 
11672.81 



0.02429 
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IV. COMPARISON OF ANALYTICAL AND COMPUTER SIMULATION RESULTS 



By comparing the analytical results with the simulation 
results, it is possible to evaluate the analytical model. 

If the measures of effectiveness predicted by the analytical 
model are reasonably close to those generated by the simula- 
tion, there is support for the analytical results. In the 
table below we compare the main measures of effectiveness: 
back orders and costs. 

Total Back Orders Total Variable 
Per Year Costs Per Year 



Analytical result 175.56 $10,808.17 

Simulation result 390.00 $11,672.81 or 

(For both starting 

conditions) 256.1 $11,389.25 

This table shows that the average yearly back orders are 
considerably higher than what is estimated by the mathematical 
solution. In fact, the total number of back orders generated 
in the simulation does not even meet the constraint. Since 
the solutions produced by the simulation model, especially 
with respect to the number of back orders do not agree with 
the mathematical model, this suggests that some of the 
assumptions made in deriving the mathematical results are not 
reasonable. The major differences between the simulation 
results and the analytical results are found in the measures 
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of effectiveness at the main system. The results for the 
lower echelon track reasonably closely. 

The costs are also a little higher than what is 
expected but much closer percentage wise. 

As an explanation for the large differences observed 
in the number of back orders at the main system with the 
simulation and analytical models, let us reexamine the 
assumptions we made that affect back orders. Back orders 
results from inadequate amounts of safety stock to protect 
the inventory system against excessively large numbers of 
demands in a lead time. The safety stock is manipulated 
by control of the reorder point r. After the reorder point 
is hit and an order placed, the system is totally at the 
mercy of the demands that occur during the lead time. If 
the variance of lead time demand is underestimated, then 
large stockouts will occur, even if the mean lead time 
demand is estimated accurately. The lead time demand is 
affected by two random quantities: (1) the demand distri- 

bution and (2) the lead time distribution. 

For purposes of making the mathematics tractable, we 
assumed in our analytical model that demands at the main 
system flowed in at a smooth continuous rate A which was 
taken to be the sum of the rates of demand incurred directly 
at the main system A m and those which occurred at the lower 
echelon systems A^ and A 2 « The lead times were assumed to 
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be normally distributed with mean equal to Ax and variance 
= \(1-Ag) where x is the mean value of the lead time. 

Because of the relatively high demand rates used in the 
example runs, the normal assumption should be justified by 
the Central Limit Theorem. However, as is illustrated by 
the Versatec plot output (See Appendix F) the demands at the 
main system are far from smooth. What happens is the demands 
which arrive directly at the main system cause the inventory 
position to drop off smoothly. However, when the replenish- 
ment orders from the lower echelon are received, large drops 
occur in the inventory position of the main system. Recall 
that the lower echelon systems order in batches once each 
period from the main system. If the main system simply tries 
to average out demands (as assumed by the analytical model) , 
it will sometimes have very large amounts of excess stocks 
when shipments arrive and sometimes very large numbers of 
back orders. The high variance in lead time demand caused 
by the irregular demand actually seen at the main system 
causes the problem. Since all demands eventually flow 
through the main system and this is assumed by our analytical 
model, the problem cannot lie in the value used for the 
mean lead time demand. 

Let us explore further what happens when reorders are 
triggered at the main system. Demands directly at the main 
system eat away smoothly at the inventory position. Then a 
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very large quantity, say X, is demanded by one of the lower 
echelon systems. There is a very good chance that the 
large order placed by the lower echelon system will trigger 
a reorder by the main system. However, if the demand causes 
a large overshoot of the reorder point, the main system may 
have much less stock to live off of until the shipment 
arrives. For example, suppose the inventory position is 
IP = 347, the reorder point is r = 300, and System One places 
it's resupply order for 200 units. An order will be placed 
by the main system but instead of having 300 units to keep 
it going until the order is received, it will have only 
147 units. It is clear that if 300 is the amount of stock 
needed to provide reasonable protection against demands in 
a lead time, large numbers of back orders would be expected. 
Effectively, in the example above, the reorder level was not 
r = 300 but r = 147 and the safety stock negative. The 
impact of this surge in demand caused by the batching of 
demands received from the lower echelon systems is to reduce 
the "effective" reorder point from the value r to a value 
r'<r. Because the actual demand distribution witnessed b, 
the main system is difficult to describe mathematically, 
the actual value of r' cannot be determined analytically. 
However, it is clearly less than r and may be much more so. 

In the next chapter we describe a policy modification 
for operation of the main system that was suggested by the 
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observations above. The modification attempts to allow 
the main system to anticipate the surge of demands that 
will be received by the lower echelon systems. 
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V. AN ALTERNATIVE SOLUTION 



A. DESCRIPTION OF ALTERNATIVE PROCEDURE 

Due to the reasons mentioned in Chapter IV, the number 
of back orders found by simulation were much higher than 
the number of back orders predicted by the mathematical 
model. The question is how to run the multi-echelon system 
so that the large number of back orders seen earlier can 
be reduced. 

Obviously, what needs to be done is to reduce the impact 
of the very large demands that occur when the lower echelon 
systems place their resupply orders. With modern day 
communication and data systems, it would be feasible to 
allow the main system to "see" every demand that occurs 
anywhere in the system. If the main system is given the 
visibility, it will be able to take action to get the stock 
on its shelves in anticipation of the large demands from 
the lower echelon systems. The main system can do this if 
it uses as its reorder point the pseudo inventory position 
which is like the inventory position except that it decreases 
only when direct customer demands are encountered at any 
of the three systems. The pseudo inventory position is un- 
affected by the batch replenishment demands placed by the 
lower echelon systems. For example, if a customer requests 
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a unit from System One, the pseudo inventory position at the 
main system decreases by one. However, if System One places 
an order for 200 units at the main system, the pseudo inven- 
tory position does not change. 

This new policy is referred to in this thesis as the 
"early warning policy." Clearly, the pseudo inventory posi- 
tion will always reach the reorder point before the inventory 
position. Therefore, orders will always be placed earlier 
and consequently, the number of back orders should decrease. 
The price paid will be in terms of extra holding costs. The 
results of a simulation using the early warning policy should 
be more nearly like those of the mathematical model since, 
in effect, the mathematical model makes the assumption of 
early warning. Implicitly, assumption of demand which is 
the superposition of the direct demands occurring at the 
main system and Systems One and Two is equivalent to the 
"early warning assumptions." 

In the next section, results are given of a simulation 
of the multi-echelon system with early warning. The flow- 
chart of the simulation model is given in Appendix B. The 
actual FORTRAN computer program is given in Appendix E. 

B. COMPUTER SIMULATION RESULTS FOR EARLY WARNING POLICY 

The same starting conditions and input parameters used 
in the previous simulation run were utilized here. The 
results are summarized as follows: 
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1 . 



Main System Results 



The n limber of orders in 10 years = 50. 

Total number of back orders in 10 years = 126. 

The average on-hand inventory = 488.093. 

Average safety stock = 210.96. 

The costs are: 

Ordering costs = $2,500.00 
Holding costs = 5,613.00 

Back order costs = 63.00 

Total costs $8,176.07 

2 , System One Results 

Total number of back orders is 10 years = 665. 

The average on-hand inventory over 10 years = 117.846. 
The costs are: 

Ordering and reviewing costs = $1,006.18 
Holding Costs = 1,355.23 

Backorder costs = 332.50 

Total costs $2,691.91 

3. System Two Results 

Total number of back orders in 10 years = 312. 

The average on-hand inventory over 10 years = 58.552. 
The costs are: 

Ordering and reviewing costs = $ 589.68 

Holding costs = 673.35 

Back orders costs = 156.00 

Total costs $1,419.03 

C. COMPARISON OF EARLY WARNING POLICY SIMULATION RESULTS 
WITH ANALYTICAL AND FIRST SIMULATION RESULTS 

The comparison will be done with respect to total 
yearly back orders and total variable system costs (con- 
sidering all three system entities) . 
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Entire System 



Number of Back Total 

Orders Per Year Cost 

Analytical Result 175.56 . $10,808.47 

First Simulation Results 256.1 11,389.25 

Early Warning Simulation 

Results 110.3 12,287.01 

Since the effects of implementing the early warning policy 
will be observed primarily at the main system, we also 
produce the results obtained for the main system individually. 



Analytical Result 

First Simulation Result 

Early Warning Simulation 
Results 



Main System 


Numbers of Back 
Orders Per Year 


Total 

Cost 


108.05 


$6,796.50 


166.50 


7,322.28 


12.6 


8,176.07 



This table shows the differences better than the first one. 
The number of back orders decreases 88.3 percent, simultaneously 
as the costs increase about 10.44 percent. The reason for the 
higher cost is because the increment in safety stock increases 
the carrying costs. 
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VI. CONCLUSIONS 



The differences between the results of the mathematical 
and the simulation model can be explained largely as a result 
of the assumptions made about the demand process. In the 
main system it was assumed that the demand was smooth (the 
superposition of three Poisson processes) , but in the simu- 
lation model the actual demand at the main system was as it 
would be in actual practice. There were batches of demand 
placed by System One and System Two in addition to the indi- 
vidual customer demands directly at the main system. These 
demand batches in fact increased the variance in lead time 
demand beyond that modelled. This explains why many more 
back orders were generated in the simulation model than 
what was predicted by the mathematical model. 

The simulation model developed in this thesis is useful 
for making comparisons and examining the effects of policy 
changes or parameter changes. Moreover, it is one of the 
best ways to check the reasonableness of all of the simpli- 
fying assumptions made in order to obtain analytical solutions. 

The mathematical results described in this thesis do 
not adequately determine the reorder point. The predicted 
number of stockouts is much less than the simulated numbers. 

As explained earlier, this is probably due to the assumptions 
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made in the model about the variance of lead time demand. 
Further study is needed to determine what could be done 
in the analytical model to better approximate what happens 
in actual practice. 

The early warning policy discussed in this thesis did 
provide for great reductions in the number of stockouts 
system wide. Since stockouts are probably the most important 
consideration in military supply systems, the early warning 
policy is recommended, even though the holding costs are 
larger. Additional study is required to see if the reorder 
level in the early warning policy can be reduced substantially 
from the value determined by the mathematical model. Pre- 
liminary evidence is that the reorder level can be reduced 
significantly (25 percent or so) without generating excessively 
many back orders if the early warning policy is used. 

Our simulation models allow us to view the effect of 
changes but they cannot be used to optimize the values of 
the policy parameters. For that objective, additional work 
in the mathematical modelling area is required. 

In this thesis we have tried to model a multi-echelon 
inventory system analytically, by linking together the 
individual echelons and locations through a single objective 
function and a constraint on backorders system wide. We, 
knowingly, were making various simplifying assumptions to 
facilitate the derivation of solutions. As reported above. 
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the resulting solution for the reorder level at the main 
system led to many more back orders than predicted. The 
other solutions; , R2, Q, T^, and T2 appear to be satis- 
factory. 

In order to model adequately the multi-echelon system, 
it will probably be necessary to build into the determina- 
tion of the reorder lead at the main system the values of 
the parameters R^ and R2 at the lower echelon systems. 

This will be the only way to accurately describe the actual 
demand process that is observed at the main system. We 
recommend future work in this area. 
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APPENDIX A 



FLOWCHARTS OF FIRST SIMULATION MODEL 
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Subroutine 

Subroutine 

Subroutine 

Subroutine 

Subroutine 

Subroutine 



One = Ship arrivals to System One 
Two = Ship arrivals to System Two 
Three = Ship arrivals to Main System 
Four = Periodic review of System One 
Five = Periodic review of System Two 
Six = Shipment arrival to System One 
Seven = Shipment arrival to System Two 
Eight = Shipment arrival to Main System 
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APPENDIX B 



Flowcharts of Early Warning Simulation Model 
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Subroutine 
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One = Ship arrivals to System One 

Two = Ship arrivals to System Two 

Three = Ship arrivals to the Main System 

Four = Periodic review of System One 

Five = Periodic review of System Two 

Six = Shipment arrival to System One 

Seven = Shipment arrival to System Two 

Eight = Shipment arrival to the Main System 
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APPENDIX C 



Simulation Program; Versatec Plotter 
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APPENDIX F 

Versatec Output of the Program in Appendix C 
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APPENDIX G 



TI-59 Calculator Program for Analytical Solutions for 
Periodic Review Systems 

This program was written for probabilistic periodic 
review inventory models having gamma distributed lead times 
and Poisson arrivals and uses. Normal distribution for 
demand during lead time instead of Negative Binomial. 

Input Requirements 

The following variables should be stored in the regis- 
ters shown before the variables. 

STO 01 = A (arrival rate per day) 

STO 02 = C 
STO 03 = I 
STO 04 = n 
STO 05 = J 
STO 06 = A 
STO 07 = a-1 
STO 08 = 6 

This program also requires the TI-59 applied statistics 
module . 
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APPENDIX H 



Variable Definitions for Simulation Programs 



IQ 


= Net inventory 


IP 


= Inventory position 


IPP 


= Second inventory position for early warning system 


IR 


= Order quantity 


IS 


= Reorder level 


ID 


= Number of items demanded by systems 


IF 


= Number of items demanded per demand by ships 


QU 


= Ship queue 


QS 


= Lower echelon's group demand queue in the main 
system 


IG 


: Amount of demand for each demand waiting in the 
queue 


IA 


= Amount of shipment arrived 


IM 


= Index. If it is equal to 1, that means that the 
demand waiting in the main system queue to be 
filled belongs to System One. 


TW 


= Total waiting time 


TOH 


= Total average on-hand inventory 


T 


= Length of a period 


EVENT = 


= This indicates the subroutines 


I EVENT = 


= This indicates time of subroutines scheduled 


X 


= Net inventory variable for Versatec plotter 


V 


= Inventory position variable for Versatec plotter 


Y 


= Time for Versatec plotter 
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WK = Work space for geometric random variable 

WS = Work space for geometric random variable 

WZ = Work space for geometric random variable 

S = Exponential random number 

IK = Geometric random number 

WW = Increment 

I = Indicates the number of ships waiting in the ship 

queue at System One 

J = Indicates the number of ships waiting in the ship 

queue at System Two 

L = Indicates the number of ships waiting in the ship 

queue at the main system 

K = Indicates the number of demand batches waiting 

in the group demand queue at the main system 

SS = Time indicator 

IQQ = Counter for ship arrivals to indicate the number 
of items demanded per period 

ML = Multiplier for the number of batches of demand to 

be ordered from outside supplier 

KA = Indicates the last change on net inventory of the 

main system in order to get the number of items 
demanded in a lead time 

AK = Time of last change 

A = Scale parameter for a gamma distribution 

B = Shape parameter for a gamma distribution 

XM = Ship arrival rate per day 

P = Probability of success for geometric distribution 
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